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Shifts  of  the  intra-ligand  stretch  frequencies  of  CO  and  CN  chemisorbed  on  metal 

1-3 

surfac  .  lave  been  observed  when  an  electric  field,  F,  is  applied.  Experiments  have 

been  carried  out  for  CO/Ni(110)  in  UHV,1  for  CO  on  Pt  electrodes  in  an  electrochemical 
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cell,  and  for  CN  on  several  metals  in  v  .ctrochemical  environments.  Two  distinct 

explanations  for  the  field  induced  vibrational  shifts  have  been  proposed.  Lambert1  has 

argued  that  the  shift  is  a  Stark  effect  with  little  enough  changes  in  the  metal-ligand 

chemical  bonding  that  it  can  be  described  by  perturbation  theory.  On  the  other  hand, 

4 

Holloway  and  Norskov  have  argued  that  the  field  changes  the  nature  of  the  chemical 

bond  and  the  frequency  shift  results  from  changes  in  the  CO  2w  occupation.  In  this 

letter,  we  present  cluster  model  studies  of  CO/Cu(100)  and  CN/Cu(100).  For  the 

changes  in  vibrational  frequency,  u,  due  to  F,  we  find  that  the  Stark  perturbation  theory 

treatment  describes  the  shift  of  the  u  rather  well  and  that  chemical  changes  are  small. 

Further,  the  most  important  effect  of  the  field  is  to  change  the  equilibrium  bond  length, 

re,  which  couples  with  the  anharmonicity  of  the  potential  energy  curve  to  cause  a  change 

in  u.  However,  since  the  metal-ligand  bonding  is  quite  different  for  CO  and  CN,  the  way 

in  which  the  field  affects  «  is  also  different.  For  CO,  the  covalent  bonding,  dominantly 
*  5,6 

metal  to  CO  2v  back-donation,  leads  to  a  direct  change  in  the  ligand  a.  For  CN  where 
the  bonding  is  mainly  ionic,  the  direct  change  in  the  ligand  u>  is  rather  small  while  the 
change  in  the  metal-CN  u  is  large.  The  coupling  of  these  stretching  modes  leads  to  the 
change  in  the  higher  frequency,  dominantly  C-N  stretch,  normal  mode. 

We  describe  first  the  clusters,  adsorption  geometry,  electronic  wave  functions,  and 
internal  coordinates  used.  Next,  we  consider  the  nature  of  the  metal-ligand  bond  in  the 
absence  of  a  field.  Finally,  we  show  how  the  ligand  u  is  changed  when  an  electric  field  is 


applied. 


2 


The  cluster  used  to  model  Cu(100)  has  10  atoms,  five  in  the  first  layer,  four  in  the 
second  and  one  in  the  third;  the  bond  distances  are  from  bulk  Cu.  The  ligand  is  added  at 
an  on  top  site  with  C  nearest  the  surface  and  the  ligand  axis  normal  to  it.  For 
C0/Cu(100),  as  well  as  other  metals,  there  is  strong  evidence  for  this  adsorption 
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geometry.  For  CN/Cu(100),  this  linear  geometry  is  often  found  for  metal-CN 
8 

complexes.  This  cluster  and  others,  larger  and  smaller,  have  been  used  for  the  study  of 
6.9,10 

several  properties  of  CO/Cu(100).  Self-Consistent  Field,  SCF,  cluster  wave 

functions  are  obtained  for  several  Cu-C  and  C-0  (C-N)  bond  distances,  with  and  without 

of  a  uniform  electric  field  normal  to  the  surface.  The  SCF  wave  functions  use  a  pseudo 

potential  for  the  core  electrons  of  the  nine  environmental  Cu  atoms.  Refs.  10  and  1 1 

discuss  the  computational  approach  in  more  detail.  For  Cu1QCO,  the  ligand  geometry  is 

varied  about  the  positions  R(Cu-C)=3.70  and  R(C-0)=2.15  bohrs:  these  are  close  to  the 

equilibrium  distance  from  low-energy  electron  diffraction7  and  from  calculations  for  a 
12 

Cu-CO  cluster.  For  Cu1()CN,  the  equilibrium  R(Cu-C)  and  R(C-N)  were  determined 
and  the  geometry  varied  about  bond  distances  near  equilibrium.  For  the  ligand  vibration, 
an  internal  coordinate  which  fixed  the  ligand  center  of  mass  was  used.  For  the  metal-CN 
stretch,  the  CN  was  translated  with  fixed  R(C-N).  The  equilibrium  bond  lengths,  r  ,  and 
the  u  were  determined  from  a  polynomial  fit  to  the  energy  curves  along  these  internal 
coordinates.  Even  though  our  computational  approach  cannot  be  expected  to  yield  exact 
absolute  values  of  w,  the  relative  changes  due  to  the  electric  field  will  be  well  represented. 

Our  results  show  that  the  bonding  of  CN/Cu(100)  is  dramatically  different  from 

CO/Cu(100).  We  use  the  projection  of  the  free  ligand  orbitals,  «t>p,  on  the  cluster 

wavefunction  to  characterize  the  transfer  of  charge  between  the  metal  and  ligand  units. 

This  projection,  which  gives  a  more  realistic  view  of  the  charge  transfer  than  a 
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population  analysis,  is  summarized  in  Table  I. 


The  charge  transfer  between  a  metal  and  CO  is  normally  divided  into  a  donation 
from  CO  to  the  metal  and  back-donation  from  the  metal  to  the  unoccupied  CO  2n.  The 
projection  of  the  CO  1  a  to  5 a  and  lir  orbitals  measures  the  extent  to  which  these  orbitals 
are  unoccupied  in  the  cluster.  For  Cu^CO,  this  projection  gives  nearly  14  electrons  in 

these  CO  orbitals,  see  Table  I,  showing  that  the  CO  donation  is  small.  The  projection  of 

•  * 

the  CO  2 n  orbital  measures  the  metal  back-donation.  Approximately  0.2  2w  electrons 

are  occupied  in  Cu10CO  bringing  the  total  number  of  electrons  on  CO  to  14.2.  While  this 
is  not  a  precise  measure  of  the  CO  ionicity,  it  does  indicate  the  build-up  of  it  charge 
between  CO  and  the  Cu  surface  expected  for  a  covalent  bond. 

The  projections  for  the  [Cu1()CN]  ,  [Cu^CN]0,  and  [CuJQCN]+  clusters  confirm 
that  CN  is  adsorbed  essentially  as  CN  .  In  all  three  cases,  the  projection  of  the  occupied 
orbitals  of  neutral  CN  indicates,  as  for  CO,  that  almost  14  electrons  occupy  these 
orbitals.  A  further  test  of  the  CN  character  is  made  by  projecting  the  free  CN  orbitals 
from  the  Cu1()CN  wave  function;  this  leads.  Table  I,  to  a  projection  even  closer  to  14.  The 
projection  of  CN  2a-  from  Cu^CN  is  smaller  by  one-half  than  for  Cu10CO  indicating 
a  smaller  v  back-donation  to  CN.  We  consider  only  the  Cu)()CN  stretching  frequencies 
in  the  following. 

The  field  free  CO  bonding  involves  an  almost  neutral  ligand  with  a  significant  dative 
covalent  bond  between  the  metal  and  the  ligand  2ir;  the  CO  a  donation  is  very  small.  In 
contrast,  the  CN  bonding  on  a  Cu  surface  is  largely  ionic;  the  covalent  bonding,  a 
donation  and  tt  back-donation,  is  small. 

For  a  cluster  in  a  uniform  electric  field,  F,  the  SCF  first  order  perturbation  theory 
energy,  Ep(F),  is  obtained  by  subtracting  n  •  F  from  the  field  free  SCF  energy  where  g 
is  the  field  free  cluster  dipole  moment.  The  SCF  variational  energy,  E„r  (F),  is  obtained 


by  adding  (Sr^  -  SQ^-R^)  •  F  to  the  field  free  Hamiltonian  where  i  denotes  electrons  and 


N  nuclei;  Q..  is  the  nuclear  charge.  These  E^fF)  and  EC__(F)  are  for  given  nuclear 

positions.  The  fields  considered  are  F=  ±0.01  a.u.  =  ±5.2  x  107  V/cm.  For  reference, 
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the  field  used  by  Lambert  in  his  UHV  study  of  CO/Ni(110)  was  6.0  x  10  V/cm.  In  an 

electrochemical  cell,  an  externally  applied  potential  of  1  V  can  create  a  field  gradient  of 

~107  V/cm  at  the  electrode  surface.15  Thus,  a  field  of  F=0.01  a.u.  though  considerably 

larger  than  used  by  Lambert1  is  comparable  to  those  believed  to  exist  at  electrochemical 

interfaces.  Values  of  w  and  the  equilibrium  bond  length,  rg,  computed  from  the  Ep(F)  and 

Escf(F)  curves  along  an  internal  coordinate  are  given  in  Table  II. 


For  the  CO  stretch,  the  Stark  and  SCF  frequency  shifts,  Aw,  are  the  same  within 

5%.  Regarding  F  to  be  a  variable,  the  average  tuning  rate,  dw/dF  =  1.3  x  lO"6 

cm  V(V/cm),  is  close  to  the  1.1  ±0.4  x  10"6  measured1  for  CO/Ni(110).  For  the  CN 

stretch,  the  Aw  for  both  Stark  and  full  variational  w’s  are  very  small.  This  is  interesting 

3  2 

since  chemisorbed  CN  has  a  similar  measured  tuning  rate  as  CO.  However,  the  tuning 
rate  for  the  metal-CN  stretch  is  quite  large.  It  is  still  largely  a  Stark  effect.  We  show 
below  that  the  coupling  of  the  C-N  and  Cu-CN  internal  modes  leads  to  a  larger  shift  in 
the  high  frequency,  dominantly  C-N,  normal  mode. 

Lambert1  made  Taylor  series  expansions  of  the  molecular  potential  and  the  dipole 
moment  curves  to  obtain  the  first  order  Stark  shift.  We  carry  out  this  expansion 
somewhat  differently  to  stress  the  two  different  effects  involving  the  field.  The 
perturbation  theory  potential  curve  is: 


V(r,F)  -  ao  -  MqF  -  M,F(r  -  re) 1  +  (a2  -  M2F)(r  -  re)2  +  .... 

■ 


* 

1“ 

r 


(i) 


where  the  expansion  is  about  the  zero  field  rg  and  the  expansion  coefficients  for  the 
potential,  a.,  and  the  dipole  moment,  M  ,  are  for  the  F=0  wave  function.  The  first 
derivative  of  the  dipole  moment,  Mj,  leads  to  a  change  in  rg: 

Are(F)  =  M,F/2(a2  -  M2F)  .  (2) 

The  change  in  the  curvature,  d2V(r,F)/d2(r  -  re),  at  r  (F)  has  a  direct  contribution  from 
the  curvature  of  the  dipole  moment,  M2,  and  an  indirect  contribution  from  the  change  in 
r  .  The  changes  in  w,  Aw,  due  to  these  changes  in  the  curvature  are: 

(a3^e)  (3) 

Aw2  =  -^-wc(M2F/a2)  (41 

where  AWj  is  due  to  the  change  in  rg  and  Aw^  is  due  to  the  curvature  of  the  dipole  moment. 
The  field  free  frequency  is  given  by  w  and  a3  represents  the  anharmonicity  of  the 
molecular  potential  curve.  These  formulas  lead  to  a  total  Aw  =  Awj  +  Aw^  essentially  the 
same  as  Lambert’s  .  The  Taylor  series  coefficients  and  the  changes  in  re  and  w  are  given 
in  Table  01. 

For  the  ligand  stretches,  the  exact  perturbation  theory  results  of  Table  II  are  closely 

reproduced.  The  value  of  Aw.,  is  much  larger  than  that  of  Aw„.  The  Ar  and  Aw,  are  much 

1  2  e  1 

larger  for  the  CO  stretch  than  for  the  CN  stretch  because  the  first  derivative  of  the  dipole 
moment  curve,  M]t  is  larger,  by  a  factor  of  6,  for  CO.  For  free  CO,  M^-1.02.  The 
Cu10CO,  or  C0/Cu(100),  Mj  is  almost  twice  as  large  as  the  free  CO  value.  We  have 
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determined  that  the  Cu  to  CO  w  back-donation  makes  the  largest  contribution  to  this 
increase. 

For  free  CN,  the  small  Mj=  -0.4  a.u.  arises  because  it  is  isoelectronic  to 
homopolar  N2  which  has  ju.  =  0.  (For  CN  ,  g  is  taken  for  a  point  fixed  with  respect  to  the 
CN  center  of  mass.)  The  Cu-CN  bond  is  largely  ionic  but  the  orbital  projection.  Table  I, 
indicates  small  a  donation  and  it  back-donation  contribute  to  the  decrease  of  |  M.  |  for 
the  Cu1QCN  where  Mj  =  -0.3  a.u..  A  larger  cluster  better  able  to  represent  the  image 
charge  could  lead  to  a  larger  value  of  i  Mj  | .  The  upper  limit  is  given  by  image  charge 
theory  to  be  M^-0.8  but  this  neglects  the  chemical  changes  which  our  cluster  results 
show  reduce  j  M1 1  . 

For  the  metal-CN  stretch,  the  dipole  derivative,  is  large,  Mj  =  -1.3.  In  this  internal 
coordinate,  we  are  moving  a  negative  ion,  CN  ,  with  respect  to  the  surface.  Image  charge 
theory  predicts  a  value  of  Mj  =  -2  since  when  CN  is  translated  by  A,  its  image  also  moves 
A.  However,  even  with  the  present  cluster,  the  Taylor  series  expansion  gives  large 
|  Are  |  =  0.2  bohrs  for  |  F  |  =  0.01  a.u.  ,see  Table  HI.  These  Arg  are  smaller  than  the  exact 
perturbation  theory  results,  Table  II,  because  the  expansions  in  Eqs.  (1)  -  (4)  are  not 
appropriate  for  large  changes  in  rg  and  u.  It  is  clear  that  the  large  |  |  is  the  main 

reason  for  the  large  values  of  Au  in  Table  II. 

The  internal  C-N  and  Cu-CN  modes  are  coupled  and  the  normal  coordinates  are  a 

mixture  of  them.  In  order  to  investigate  this  coupling,  we  have  used  standard 
16 

expressions  for  the  coupling  of  two  spring-like  force  constants  which  enables  us  to 
estimate  the  magnitude  of  the  coupling.  The  differences  between  the  internal  coordinate 
frequencies  and  the  approximate  normal  mode  frequencies  are  given  in  Table  IV.  There 


is  a  clear  difference  between  the  shifts  in  the  higher  frequency  modes  for  different  values 


of  F.  When  this  effect  is  added  to  the  Aw  in  Table  II,  we  obtain  Aw(F  =  +0.01)  =  15  and 
Aw(F  =  —0.01)  =  —24cm-1  significantly  larger  values  than  those  for  the  C-N  internal 
coordinate.  Thus,  the  shift  of  the  high  frequency  normal  mode  is  because  there  is  a  shift 
in  the  metal-CN  stretch  which  is  coupled  with  the  ligand  stretch.  However,  our  present 
treatment,  cluster  size  and  coupling  of  the  internal  modes,  cannot  give  quantitative  values 
for  the  shift. 

For  Cu^CO,  we  have  calculated  the  dipole  moment  derivative,  Mj,  near  the 
metal-ligand  separation,  R(Cu  —  C)  =  3.70  bohrs,  used.  We  find  |  |  <  0.05.  For  the 

low  frequency,  metal-CO,  normal  mode  for  CO/Ni(100),  the  value  Mj  =0.2  has  been 
measured.17  This  small  value  means  that  the  metal  ligand  internal  coordinate  stretch 
frequency  will  hardly  change  in  an  electric  field.  Accordingly,  the  shifts  due  to  the 
coupling  of  the  internal  coordinates  will  be  nearly  constant  irrespective  of  imposed 
electric  field. 

In  conclusion,  we  have  shown  that  the  effect  of  an  electric  field  on  the  stretching 
frequencies  of  CO  and  CN  chemisorbed  on  Cu(100)  is  described,  to  a  very  good 
approximation,  by  perturbation  theory.  The  most  important  effect  of  the  field  is  to 
change  the  equilibrium  bond  length.  This  change  coupled  with  the  anharmonicity  of  the 
potential  is  the  dominant  reason  for  changes  in  the  w.  Because  the  metal-ligand  bonding 
is  very  different  for  CO  and  CN,  the  effect  of  the  field  is  different  is  these  two  cases.  For 
CO,  the  covalent  chemical  bond  leads  to  large  changes  in  w  for  the  C-O  stretch  internal 
coordinate.  For  CN,  the  bonding  is  ionic,  chemisorbed  CN  is  essentially  CN  .  The 
derivative  of  the  dipole  moment  for  the  CN  stretch  internal  coordinate  is  small  and, 
hence,  the  change  in  w  is  small.  On  the  other  hand,  the  dipole  derivative  for  the  metal-CN 
stretch  is  large  and  this  is  transferred  to  the  CN  stretch  through  the  coupling  of  these  two 


coordinates.  We  predict  that  there  should  be  a  large  field  induced  shift  in  the  low 

frequency  metal-CN  stretch  while  the  shift  of  the  metal-CO  stretch  should  be  small.  This 

could  be  tested  by  Surface  Enhanced  Raman  Spectroscopy  in  which  the  low  energy  modes 
3b 

are  readily  accessible  .  These  origins  of  the  field  induced  a  shifts  are  for  Cu  however 
we  expect  them  to  be  true  for  other  metal  surfaces  too,  in  particular  the  noble  metals, 
since  the  ligand  bonding  is  likely  to  have  similar  common  features. 

This  work  has  been  supported  in  pan  by  the  Office  of  Naval  Research. 
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TABLE  I 


Projections  of  the  occupied  and  2 ir*  ligand  orbitals  on  the  Cu1QCO  and  Cu^CN  clusters.  For  the  occupied 
orhitals.  the  projections  from  the  a  and  r  orbitals  are  given  separately.  P0  and  P_.  For  Cu  CO. 

R(Cu-C)  =  3.70  and  R(C-0)  =  2.15  bohrs:  for  CuJQC\.  R(Cu-C)=3.90  and  R(C-N)  =  2.20  bohrs. 


Ligand  Orbitals 

Projections 

Cu10 

CO 

[Cu10CNjO 

O 

u 

CN  ]  + 

O 

3 

U 

CN]- 

L;  occupied 

P, 

9.99 

P, 

4.00 

P,  P» 

9.92  3.98 

p„ 

9.91 

u> 
vO  ri 
oc 

p„ 

9.91 

P* 

3.98 

L';occupiea 

— 

— 

9.96  3.99 

9.96 

3.98 

9.97 

3.99 

L";2t" 

— 

0.23 

0.14 

— 

0.13 

— 

0.14 

11 


TABLE  III 

Taylor  series  expansion  coefficients,  in  a.u.  for  the  molecular  potential,  a.,  and  dipole 
moment.  \I ,  curves  along  the  ligand  and  metal-ligand  coordinates.  The  perturbation  theory 
changes  for  r  ,  Ar  in  bohrs,  and  for  frequency.  Aw  ,  Aw  ,  and  Aw  =  Awi  +  Aw2  in  cm  1  are 
also  given. 

Stretch  C-O  C-N  Cu-CN 


TABLE  IV 


Differences,  Au,  between  the  approximate  normal  modes  and  the  internal  coordinates  for 
CUioCN. 


F(a.u.) 


Aojfcm'1) 


Lower 

Mode 


Higher 

Mode 

+  15.6 
+26.1 
+  4.7 


m**'- 


S 
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